Holon Pairing Instability based on the Bethe-Salpeter Equation obtained from the t-J 
Hamiltonians of both U(l) and SU(2) Slave-boson Symmetries 
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We investigate a possibility of holon pairing for bose condensation based on the Bethe-Salpeter 
equation obtained from the use of the t-J Hamiltonians of both the U(l) and SU(2) slave-boson 
symmetries. It is shown that the vertex function contributed from ladder diagram series involving 
holon-holon scattering channel in the Bethe-Salpeter equation leads to a singular behavior at a 
critical temperature at each hole doping concentration, showing the instability of the normal state 
against holon pairing. We find that the holon pairing instability occurs only in a limited range of 
hole doping, by showing an "arch" shaped bose condensation line in agreement with observation for 
high Tc cuprates. It is revealed that this is in agreement with a functional integral approach of the 
slave- boson theories. 



Since the advent of high Tc superconductivity, both 
the U(l) and SU(2) slave-boson approaches to the t-J 
Hamiltonian have been proposed to study superconduc- 
tivity for the two-dimensional systems of copper oxides 
0- Hj. Unlike the U(l) slave-boson mean field theory, 
the recently proposed SU(2) theory of Wen and Lee has 
a merit of treating the low energy phase fluctuations of 
order parameters Recently we proposed a theory of 
bose condensation by paying attention to the holon pair- 
ing channel in the U(l) slave-boson theory Q, in con- 
trast to other earlier studies concerned with the single 
holon condensation. In this paper the phase fluctuations 
of the spinon pairing and hopping order parameters are 
taken care of by the SU(2) theory. In addition, follow- 
ing our earlier study the spinon and holon degrees 
of freedom are introduced into the Heisenberg term on 
the basis of on-site charge fluctuations which arise as a 
result of site-to-site electron(and thus holon) hopping in 
two-dimensional quantum systems of hole-doped high 
cuprates. Most of the slave-boson theories have been 
concerned with single holon bose condensation. These 
theories predicted a linear increase of bose condensation 
temperature with doping concentration, contrary to the 
observed bose condensation line of an 'arch' shape which 
manifests the presence of optimal doping rate. Most re- 
cently Wen and Lee |3| questioned whether the single 
holon or the holon pair bose condensation is favored with 
the SU(2) slave-boson theory. In this work we pay atten- 
tion to holon pair bose condensation by demonstrating 
a possibility of holon pairing instability from the study 
of the Bethe-Salpeter equation obtained from the use of 
the t-J Hamiltonians of both the U(l) and SU(2) slave- 
boson symmetries. The present approach is shown to 
predict the observed characteristics of the bose conden- 
sation line of an arch shape in excellent agreement with a 
functional integral approach to the slave-boson theories. 

We write the t-J Hamiltonian, 
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H = -t (cLcjT + c.c.) + JY1 ■ ^3 ^ f '^^"i)' (1) 

<i,j> <i,i> 



where S, -S,- 



Here Si is the electron spin operator at site i, = 
\c\^iy apCifj with cTq.^, the Pauli spin matrix element and 

rii, the electron number operator at site i, rii — cl^Cia- 
In the slave-boson representation 0- Q Ca, the electron 
annihilation operator of spin a can be written as a com- 
posite of spinon and holon operators. That is, = fa 
in the U(l) theory and Cq = -^h'^ipa in the SU(2) theory 
with a — 1,2, where fa{b) is the spinon(holon) annihila- 

/i 



tion operator in the U(l) theory, and t/'i 

/2 



and 



are respectively the dou- 



blcts of spinon and holon annihilation operators in the 
SU(2) theory. 

After Hubbard Stratonovich transformations in associ- 
ation with the direct, exchange and pairing channels and 
the saddle point approximation, the effective Hamilto- 
nian is decomposed into the spinon sector and the holon 
sector separately. Here to explore the instability of the 
normal state against holon pairing, we pay attention to 
the holon sector of the Hamiltonian. The holon Hamil- 
tonian is derived to be, for the U(l) slave-boson theory 



t-J,U(l) 



't XtMbj +C.C. 



<i,3> 



J 



Y H\'blb%b,-^,oJ2bh, (2) 



<jj> 



where Xtj =< fLf^o 



> is the hopping or- 
der parameter and A{j —< fjifn — fjifii >, the spinon 
pairing order parameter, and for the SU(2) slave-boson 
theory. 
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<«J> 



<ij>,Q,/3 



where J/i 



-A 



'■^ is the order parameter 

matrix of hopping order, Xij ^-nd spinon pairing or- 
der, Af^. with X»j f]J^a + Jirh)^{b]ihl - 4^J2) > 

and A{j =< fjifa - fj2fzi >■ 

Introducing a uniform hopping order parameter, Xji = 



X, and a d-wave spinon pairing order parameters, Aj^ = 
±A/(the sign +(— ) is for the ij link paraUel to x (y)), we 
obtain from Eq.(H) the momentum space representation 
of the Hamiltonian, for the U(l) theory. 



<xt/(i)=E<k)6t6k 

k 

2PJ Z-^ ^ ' -k +q k K+q, 



(4) 



k,k',q 



where e(k) — — 2tx7k — /^Oi the energy dispersion of holon 
with 7k = (cosfca; -t-cosfcy) and w(k — k) = — J|A/p7k, 
the momentum space representation of holon-holon in- 
teraction. Likewise we obtain from Eq.(^, 



ST ]J ( "*X7k - /^o tA/(/Jk \ , 
t-j,sc/(2) -2^^^\ tAf^^ txiu - MO j 

+ ^ E «(k-k>l^,^^^^&t,^6k,a6-k+q,;3, (5) 

k,k ,q,a,/3 

where a, (3 (= 1,2) represent two isospin components of 
SU(2) holon and (pk — {coskx — cos ky). N is the total 
number of sites for the two dimensional lattice of interest. 

Considering the t-matrix of involving the particle- 
particle(holon-holon) scattering channels, we obtain the 
Bethe-Salpeter equations in the U(l) slave-boson repre- 
sentation, 

<k,-k + q\t\k, -k + q >u(^i^= v(k — k) 
l-Y^ vik -k")g{k")gi-k" +q) X 



N(3 ■ 
ft 

< fc ', -fc ' -I- q\t\k, -k + q >u{i), 



(6) 



where g{k) = J^^ dre"''"^ < T^[5k(r)6[(0)] >, the holon 
Matsubara Green's function in the U(l) slave-boson rep- 
resentation, and in the SU(2) slave-boson representation, 

< k ,a ;~k + q, P \t\k, a; -k + q,P >sui2) 
^v(k' ~-k)S^, Jf^,^p 

E '"(^ ~'^')9a'a"ik")9p'p"{~k" +q) 

k" .a" ,P" 

x<k ,a ;-fc +q,(3 \t\k,a; -k + q, /3 >su{2), (7) 



where go^^ik) = dre"'"^ < T,[6ak(T)6^k(0)] the 
holon Matsubara Green's function in the SU(2) slave- 
boson representation. Here k = (fco,k) is the three 
component vector of the energy-momentum of holon. 
Q = (9o,q) is the three component vector of the total 
energy-momentum of holon pair. 

After summing over the Matsubara frequencies in 
Eqs.(^ and (j^), we obtain the following matrix equation 
for the t-matrix, in the U(l) slave-boson representation, 

E (V,k" - ™k',k" (90, q)) t^"Milo, q) = v{k - k), (8) 

k" 

where ik'.kl'JOiq) =< k',-k' +q\t\k,-k + q >u(i), 
™k',k"(9o,q) = -]^E"(^' ^^')9ik")g{-k" +q) 



i;(k' -k") 



E 



ik^ - e(k") i{-k'g + qo) - e(-k" + q) 



^ ^ " ^ .9„-(.(-k"+q)-f.(k")) ' 



and n{e) = , the boson distribution function. Simi- 
larly, we obtain, in the SU(2) slave-boson representation. 



k",a" J3" 

^v{k - k)S^-J^>p, 



(10) 



where if / '! "'^ (<Zo , q) = < k ,a ; — fc + q, (3 \t\k, a; — fc -I- 



q,P >SU{2) and m. 



a P a P' 
k',k" 



(90, q) = ^Ea;/3;^^(k 



„ 0E , (-k +q) „ 

„ n{E^,^{-k ))+e "l "(^^/(-k +q)) 

^ ) ' »go-(B„,(k") + S^^,(-k"iq)) Xf^a'a;(k Wp' f3'^{-k + 

q)^7^^„(k")[/^,^„(-k" + q). Here E,{k) = i?k - Mo 

and E2{k) — — i?k — Mo are the energy dispersions 
of the upper and lower bands of holons with i?k = 

V(X7k)2 + (A/^k)2. Uo^pik) = ~u^) 

unitary transformation matrix used for the diagonaliza- 
tion of the one-body holon Hamiltonian in Eq.(^ with 

= 75^/1-^ and = ^^1 + %^. It is noted 

that the t-matrices, ik',k(9o,q) and t^, "'^ (go , q) for the 
U(l) and SU(2) theories respectively are independent 
of the Matsubara frequencies fcg and fcp, owing to the 
consideration of instantaneous holon-holon interaction, 
w(k'-k). 
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As mentioned earlier, the t-matrix elements are ob- 
tained from Eqs.(^ and (|^, both of which involve sum- 
mation over the Matsubara frequencies. Using the ma- 
trix equations ^ and ( [lO| ) , the poles of the t-matrices are 
searched for as a function of energy qo with q = 0, i.e., the 
zero total momentum of the holon pair. The t-matrices 
are numerically evaluated from the use of Eqs.(|8|) and 
(|lO|) for each doping rate and temperature. The hopping 
and spinon pairing order parameters in Eqs.(^ and (^) 
are the saddle point values evaluated from the usual par- 
tition functions involving the functional integrals of slave- 
boson representation For both the U(l) and SU(2) 
t-matrix calculations we choose J/t — 0.3 in Eqs.(||) and 
(i- 

At high temperatures, computed poles are found to be 
positive and real, indicating that there exist no bound 
states. As temperature is lowered to a critical(onset) 
value Tc, we find that with the U(l) theory one pole 
changes its sign from positive to negative, indicating that 
their exists an instability of the normal state against 
holon pairing at the onset (critical) temperature Tc. Simi- 
larly, with the SU(2) slave-boson theory two poles change 
their signs from positive to negative at T^, these two poles 
correspond to the particle-particle(holon-holon) scatter- 
ing channels of bi-bi and 62-^2 respectively. In Fig.l. the 
onset temperature Tc (denoted by solid square in Fig.l) 
for the appearance of such negative pole(s) is predicted 
to occur only in a limited range of hole doping concentra- 
tion, by revealing an " arch" shaped feature of bose con- 
densation line in agreement with observation in the phase 
diagram of high Tc cuprates. The onset temperature 
of negative pole(s) coincide(s) surprisingly well with the 
critical temperature of holon pair condensation obtained 
from the functional integral approach(denoted by open 
square in Fig.l) of the slave-boson theories Holon- 
holon scattering occurs above the lowest possible single 
particle energy e(k = 0). The negative pole corresponds 
to the binding energy of the holon pair. This is analogous 
to the binding energy of electron pairs which results from 
the Cooper's two particle problem of electron-electron 
scattering only above the Fermi energy e(k/), i.e., the 
Fermi-surface. 

In order to find the symmetry of the holon pairing or- 
der parameter we now compute the eigenvectors of the t- 
matrix whose poles change their signs at a critical(onset) 
temperature Tc] the eigenvalue equations are 

E *k,k' (90 = 0, q = O)W^(k') = AW^(k), (11) 
k' 

for the U(l) theory, where VF(k) is the eigenvector and 
A, the eigenvalue and 

E C'''''(* = 0,q=0)I^„'^'(k') = Amk)„^, (12) 

k',Q',/3' 

for the SU(2) theory, where Wapik) is the eigenvec- 



tor concerned with the SU(2) isospin channels a and 
(3{a = 1,2 and /3 = 1, 2 ) and A, the corresponding eigen- 
value. For the case of U(l) there are N eigenvectors 
corresponding to N discrete values of momenta for an 
N X N reciprocal lattice. For SU(2), 4N eigenvectors 
are available in accordance with the isospin channels of 
holon(6Q,)-holon(&^) scattering. We search for the eigen- 
vectors whose eigenvalues diverge in the limit of go = 
and q = at an onset (critical) temperature Tc, that is, 
the eigenvectors corresponding to the pole of the t-matrix 
whose sign changes at Tc- We choose N = 10 x 10 for 
a reciprocal lattice to determine the symmetry of holon 
pairing at an underdoping rate, 5 = 0.05. 

The computed eigenvectors for both the U(l) and 
SU(2) slave-boson approaches yielded good fits to the s- 
wave symmetry of the form cos + cos ky in momentum 
space ; for the U(l) theory we have 



VF(k) = yl(cos kx + cos ky), 



(13) 



where A is the normalization constant to satisfy 
|VF(k)|2 = 1, and for the SU(2) theory, 

VF^^(k) ^ A^{coskx + cosky){Sc,i5ps + (14) 
and 

W°i^{k) = A''{coskx + cos ky) {So,, 16 PS - ^^,26^,2)^ (15) 

for bi holon and 62 holon pairing, where A*^^") is the nor- 
malization constant to satisfy X^k a /3 I^q/3°''('*)I^ ~ 1- 
One of the two computed eigenvectors shows a phase dif- 
ference of even multiples of vr between the bi and &2 holon 
pairing order parameters and is well fitted by Eq. (nj) ; the 
other computed eigenvector is well fitted by Eq.(|l5|), by 
showing a phase difference of odd multiples of tt. The 
divergence of eigenvalues was seen to occur nearly at the 
same temperature(within a numerical accuracy of lO^H). 
Fig. 2 displays the U(l) result of the s-wave symmetry of 
the holon pairing order. In Figs. 2(a) and (b), we show 
the SU(2) results of the s-wave symmetry of opposite 
phase between the 61 and 62 holon pairing order param- 
eters. Both results demonstrated excellent fits to the 
s-wave symmetry form of Eqs.(13) and (p^. Since the 
computed results which fits Eq.(14) are indistinguishable 
with the ones shown in both Fig. 2 and Figs. 3 (a) and thus 
are not displayed. In short, both the U(l) and SU(2) 
theories predicted the s-wave symmetry of holon pair- 
ing order. This is equivalent to the d-wave symmetry 
of hole pairing order, by noting that the hole is a com- 
posite of a holon(boson) of spin with charge +e and a 
spinon(fermion) of spin 1/2 with charge and realizing 
the s-wave holon pairing in the presence of the d-wave 
spinon pairing fll- |^ §. 

In the present study, we found that the bose conden- 
sation of an arch shape occurs only in a limited range 
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of hole doping concentration, by searching for the onset 
temperature of the holon pairing instabihty. It is found 
that the predicted holon pair bose condensation temper- 
ature is in precise agreements with the functional integral 
approaches of both the U(l) and the SU(2) slave-boson 
theories. Most importantly, both of these approaches 
were found to satisfactorily predict the experimentally 
observed bose condensation (superconducting) tempera- 
ture as a function of the hole doping rate, by reveahng the 
presence of the optimal doping rate in high Tc cuprates. 
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FIG. 1. Onset temperature(denoted by ■) of the negative 
energy poles representing holon pair bose condensation temper- 
ature Tc as a function of hole doping rate for both (a) U(l) and 
(b) SU(2) slave-boson theories. The symbols, Q s^nd □ represent 
respectively the spinon pairing and holon pairing temperatures 
obtained from the functional integral approach. 
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FIG. 2. The computed eigenvector of the U(l) t-matrix 
which fits the s-wave symmetry of holon pairing in Eq.(13). 
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FIG. 3. The computed eigenvector of the SU(2) t-matrix 
which fits the s-wave symmetry of holon pairing in Eq.(15). (a) 
and (b) display the phase difi'erence of tt between the 6i and 62 
holon pair order parameters. 
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